Non-perturbative vacuum-polarization effects in proton-laser collisions 
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In the collision of a high-energy proton beam and a strong laser field, merging of the laser photons 
can occur due to the polarization of vacuum. The probability of photon merging is calculated by 
accounting exactly for the laser field and presents a highly non-perturbative dependence on the 
laser intensity and frequency. It is shown that the non-perturbative vacuum-polarization effects 
can be experimentally measured by combining the next-generation of table-top petawatt lasers with 
presently available proton accelerators. 
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In classical electrodynamics the superposition principle 
holds for electromagnetic fields of arbitrary strength in 
vacuum. Instead, quantum electrodynamics (QED) pre- 
dicts that under the influence of electromagnetic fields 
of the order of the "critical" fields: E cr = m 2 c 3 /he = 
1.3 x 10 16 V/cm and B cr = m 2 c 3 /he = 4.4 x 10 13 G, 
the vacuum behaves as a nonlinear medium (— e < 
and m are the electron charge and mass, respectively) 
The experimental observation of nonlinear vacuum- 
polarization effects (VPEs) is of fundamental interest, es- 
pecially to test the validity of QED at fields of such high 
strength. Highly charged ions with a charge number Z 
such that Za < 1 with a = e 2 /Air being the fine-structure 
constant, provide electric fields of the order of E cr at the 
typical QED length scale A c = h/mc (Compton length). 
This has allowed the observation of VPEs induced by the 
strong Coulomb field of heavy ions via the measurement 
of Dclbriick scattering, i. e. the scattering of a 7-photon 
by an atomic field [2| and of the related process of 7- 
photon splitting in atomic fields [3} . Noticeably, in both 
cases effects of high-order corrections in Za have been 
observed in the cross-section 0|. 

The continuous progress in laser technology opens 
unique opportunities to test QED in the presence of in- 
tense electromagnetic waves [5(. The record laser inten- 
sity of 7 x 10 21 W/cm 2 was produced in laboratory by 
focusing a 45 TW Ti:Sapphire laser onto one wavelength 
Q. Moreover, the Optical Parametric Chirped Pulse 
Amplification (OPCPA) technique has already been suc- 
cessfully applied to generate terawatt pulses with a short 
duration (10 fs) at a high repetition rate (30 Hz) [JJ. A 
petawatt laser is currently under development based on 
the same technique Q. Finally, envisaged intensities of 
the order of 10 24 -10 26 W/cm 2 are likely to be reached 
in the near future at the Extreme Light Infrastructure 
(ELI) The first theoretical investigations of QED in 
the presence of a strong background laser field go back 



to the sixties 13, 11 1- Later, the calculations of higher- 
order QED diagrams like the mass operator , th e po- 
larization operator [13j and the photon splitting [l4| were 
performed. Refractive VPEs induced by laser fields have 
never been experimentally observed, essentially because 



present laser intensities are much smaller than the criti- 
cal intensity I cr = E 2 r /%ix = 2.3 x 10 29 W/cm 2 and laser 
photon energies much smaller than mc 2 = 0.5 MeV. To 
the best of our knowledge, the only feasible proposals to 
observe refractive VPEs in a laser field explored so far 
in the literature are restricted to perturbative VPEs in 
which the leading contributions in the laser field are the 
only, at least theoretically, observable effects while all 
higher-order effects are strongly suppressed 



15] 



In this Letter we put forward a scheme for which it is 
experimentally feasible to observe non-perturbative re- 
fractive VPEs in laser fields. We calculate the proba- 
bility for an even number of laser photons to merge in 
the head-on collision of a high-energy proton beam and a 
strong laser field due to VPEs, by accounting exactly for 
the laser field and by focusing on the non-perturbative 
dependence of this probability on the laser intensity Jo 
and frequency c^o- As a result, non-perturbative QED 
can be tested at the opening of multiphoton vacuum po- 
larization channels. 

To this end the proton beam has unique features. On 
the one hand, the proton is light enough to be acceler- 
ated up to very high energies, namely to 980 GeV at 
the Tevatron, or even to 7 TeV at the Large Hadron 
Collider (LHC) 0. Then, the laser field amplitude 
and frequency in the proton rest frame are enhanced by 
the Lorentz 7-factor of the proton so that the nonlin- 
ear QED regime becomes reachable, allowing in princi- 
ple for the experimental observation of high-order non- 
perturbative VPEs. On the other hand, the proton mass 
M = 938 MeV is large enough that the background ra- 
diation arising from multiphoton Thomson scattering of 
the laser photons is, as we will see, very low. Moreover, 
since the proton charge number is Z = 1, here the VPEs 
are essentially driven by the strong laser field. 

Without loss of generality we can consider a proton 
moving along the negative y direction with velocity (3 
and a counterpropagating laser beam with amplitude 
Eq = yJ8irIo and frequency loq (from now on natural 
units with h = c = 1 are used). The laser is assumed 
to be linearly polarized along the z direction. The pro- 
cess of laser photons merging due to VPEs in a proton 
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FIG. 1: Feynman diagram corresponding to the process of 
laser photon merging induced by VPEs in a proton field. 



field is represented by the Feynman diagram in Fig. 1 
where the right photon leg represents the outgoing pho- 
ton resulting from the merging, the thick fermion loop 
indicates that the fermion propagators are calculated 
by accounting exactly for the laser field and where the 
crossed photon leg indicates the electromagnetic field of 
the proton taken into account as a perturbation since 
Z = 1 and described by the four-potential A^(t,r) = 
eu^/4n^x 2 + j 2 (y + fit) 2 + z 2 with u* 4 = (7, 0,-/37,0) 
being the proton four- velocity and 7=1/ y/l — (3 2 its 
relativistic Lorentz factor. The probability amplitude of 
the laser-photon merging process can be calculated as 
T{k) = (2tt)- 4 / d A qA^q)T^(q,k)et(k)/V2^ (see Fig. 



1) where A^(q) is the four-dimensional Fourier transform 
of A „ (t, r), T^ u {q, k) is the polarization tensor in a laser 
field [l3[ and e v {k) is the polarization four- vector of the 
outgoing real photon with four-momentum k^ = (w,k) 
171 ]. Starting from the above amplitude, one obtains 
the differential rate d7Z by applying Fermi's golden rule: 
dK = |T(/c)| 2 dk/(27r) 3 (the quantization space- volume 
and time-interval are set equal to unity here) [13] ■ In 
the considered case of a monochromatic plane wave the 
total rate can be written as dlZ = SnLi d!Z n where dlZ n 
is the rate of photons resulting from the merging of 2n 
laser photons (Furry's theorem prevents the merging of 
an odd number of laser photons in this process [13] )■ The 
integration over k is easier by employing spherical coor- 
dinates (ui, $, ip) with the y-axis as polar axis. The inte- 
grals over u> and <p can be carried out and the ensuing 
differential rate dlZ n /d$ written as 

dK n = a 3 (l + /3)m 4 sin 3 g |ci,»| 2 + |c 2 ,„| 2 
dtf Mir 2 (1-cost?) 4 n 3 

In this expression t? is the angle between the momentum 
of the emitted photon and the y-axis and 



dv 
-1 Jo 
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;Jn(zn) ~ sin 2 pJ n {z n ) 
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where j G {1:2}, Sjji is the Kronecker S- function, 
$„ = 2np + 4p{l + e[l - sin 2 (p)/p 2 ]/2}/?7„(l - v 2 ), 
z n = 2p( 2 [sin 2 {p)/p 2 - sin(2p)/2p]/77„(l - v 2 ), A = 
1 + sm 2 (p)/p 2 - sm(2p)/p and J n {z) = J n {z) + W n {z) 
with J n {z) being the ordinary Bessel function of order 



the parameter \ r , 

Xn = £,11n : 



J 2n(l + (3)w 1 — cos $ 



(3) 



I cr m 1 + P COS 1? 
The leading-order contribution to the coefficients Cj <n de- 
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where x n — x„A A (l — v ) /96 and 
,1 
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&J> = jXn^ 



n and J' n {z) its derivative. The coefficients Cj,„ de- pends only on Xn and is given by 
pend on the two Lorentz- and gauge-invariant quantities 
£ = eEo/mcuo and rj n — (kok n )/m 2 — LUoui n (l — cos$)/m 2 
where kg = (u>o, 0, too, 0) is the four-momentum of the 
laser photons and k% is the four-momentum of the out- 
going photon whose energy w„ = 2nojo(l+j3)/(l+j3cos , d) 
depends on the number of laser photons merged, accord- 
ing to the energy conservation and the Doppler-shift. 
When the proton is highly relativistic, ui n is 47 2 times 
larger than 2nu>o for backscattered photons. 

In the asymptotic limit £ <C 1 of weak laser fields, the 
leading rate in Eq. {T]) is that with n = 1 which coincides 
with the result of the Feynman box diagram correspond- 
ing to the perturbative treatment of the laser field. The 
rate increases with raising £. Moreover, present optical 
lasers easily exceed the relativistic threshold £ ~ 1 (laser 
intensity Io w 10 18 W/cm 2 at cjq ~ 1 cV). Therefore, we 
will discuss here only the most practical regime £ 3> 1. 
In this asymptotic limit the dynamics is determined by 



f& -«p(*r/3)A-x B ^ B (4) 

A 



r7j [/ n (x„)-/;(x„)] + ^l (5) 



with /„ (x) being the modified Bessel function of order n 
and I' n {x) its derivative. As the above expressions show, 
the dependence of the differential rate d'R-n/d'd on the 
parameter \ n is highly non-perturbative. The perturba- 
tive result is obtained only if %„ -C 1 (for the process 
at hand the leading contributions at £ 3> 1 and Xn ^ 1 
coincide with those at £ <S 1 and rj n -C 1). In this case, 
the leading order asymptotics of the coefficients Cj, n have 
the typical perturbative scaling as x 2 ™: 
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2)! (2n)l3jn + l xln 



r(2n + 3/2) n! I92 r 
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with Gamma function T(x). In the particular case 
n = 1 and (3 — our expression of the total rate 
1Z\ = J dddTZ-L / ckd is in agreement with the results in 
the paper by Milstein et al. in [15|. 

In the opposite limit of large Xn we obtain 



XI 

0.136 0.201 0.316 0.525 0.834 



2e 2i7r/3 j + 1 r(2/3)r(5/6) r(n 



21 6 x /3 r(7/6) T(n + 4/3) 



1/3) 9/3 

An 
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In this regime of ultra strong laser fields where Xn 3> 1 , 
the differential rate scales non-perturbativcly with Xn, 
the higher corrections in Xn stemming from exchanges of 
laser photons with no net absorption. However, the rate 
decreases with n monotonically and behaves at large n 
as dKn/d-d ~ I/n 5 . 

We consider below some numerical examples corre- 
sponding to different possible experimental conditions in 
which the process of photon merging is feasible. In the 
first one we consider the interaction between a petawatt 
laser pulse and a proton bunch with parameters avail- 
able at the LHC. Concerning the laser, we use the fol- 
lowing data [1, @|: pulse energy 3 J, pulse duration 
5 fs at 10 Hz repetition rate. When the beam is fo- 
cused onto one wavelength of 0.8 /im its intensity is 
Io = 3x 10 22 W/cm 2 . The main parameters of the proton 
bunch are (ljj : proton energy 7 TeV, number of protons 
per bunch 11.5 x 10 10 , bunch transversal radius 16.6 /im, 
bunch length 7.55 cm. In Fig. 2 we show the differential 
rate dlZi/dd of the photons resulting from the merging 
of two laser photons (continuous line). In the upper hor- 
izontal axis in Fig. 2 we have indicated the value of the 
parameter xi at the corresponding angle it is evident 
that in the relevant region of the spectrum xx ~ 1 an d 
the perturbative approach is inadequate. For instance, 
the perturbative result for the rate dUi/dti at ■d = 3.1415 
is roughly 28 % less than the non-perturbative one. In 
this regime there is a competing process: the two-photon 
Thomson scattering of the laser photons by the proton 
beam. In fact, the energies of the photons produced via 
the photon merging and via the two-photon Thomson 
scattering are the same because of the negligible recoil 
effect in the latter case (the photon energy is 470 MeV 
at the maximal emission angle, while the proton energy 
is 7 TeV). In Fig. 2 the vacuum polarization signal 
is compared with the analogous quantity dW^'/dd of 
the two-photon Thomson scattering of the laser photons 
by the proton beam (dashed line, result via [ll|). The 
figure shows that the vacuum polarization contribution 
is clearly distinguishable from the two-photon Thomson 
scattering. By integrating dH\/dd with respect to $ and 
by multiplying the ensuing quantity with laser pulse du- 
ration and repetition rate and the effective number of 
protons that interact with the laser beam, we obtain ap- 
proximately 400 events per hour. Moreover, due to the 
non-perturbative nature of the process in this regime, 
the rate of photons resulting from the merging of jour 
laser photons is not small either: 5.4 events per hour. 
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FIG. 2: The rate per unit angle $ of photons emitted via 
VPEs (continuous line) and via two-photon Thomson scat- 
tering (dashed line). The upper horizontal axis shows the 
values of the parameter \i as a function of $ [see Eq. (3)]. 



As expected, in this case the four-photon Thomson scat- 
tering is safely negligible due to its scaling as £|, with 
£ p = (m/M)£ sa 6.4 x 10~ 2 . In the previous proposals 
|15| . the corresponding multiphoton channels involving 
more laser photons than the minimum required, are sup- 
pressed by several orders of magnitude. 

In the following example we consider the possibility 
of employing the already operative Tevatron accelerator. 
The relevant par ameters of the proton bunches at the 
Tevatron are [16(: proton energy 980 GeV, number of 
protons per bunch 24 x 10 10 , bunch transversal radius 
29 /itm, bunch length 50 cm. The relatively small pro- 
ton energy in this case can be compensated by employ- 
ing strong attosecond pulses of XUV radiation which are 
proposed in In particular, we use an attosecond 

pulse of intensity Iq = 3.8 x 10 22 W/cm 2 with a pho- 
ton energy of 70 eV and a pulse duration of 60 as. The 
latter, according to the simulations in [l8j], can be pro- 
duced by the reflection from a plasma surface of a 6 PW 
ThSapphirc laser pulse of 5 fs duration focused onto a 
spot-radius of 5 fjxa. In this setup, the angular distribu- 
tion of the photons emitted via VPEs has a shape similar 
to that in Fig. 2 and the photon rate produced via two- 
photon Thomson scattering is negligible. We report only 
the average number of photons emitted in one hour by 
the merging of two and of four laser photons, which is 
approximately equal to 490 and 6.6, respectively. As in 
[18| we have assumed that about 6 x 10 -2 of the energy 
of the initial beam is transferred to the attosecond pulse 
and that the OPCPA technique is employed to achieve 
a repetition rate of 10 Hz. Also in this case the values 
of the parameter Xn near the maximum of the photon 
angular spectrum does not permit the use of perturba- 
tion theory as xi = X2/2 « 0.7 at = 3.1410. At the 
same angle the energy of the emitted photon is about 
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440 MeV. We stress that in this example the relativistic 
parameter £ of the XUV pulse is approximately equal to 
three, which is not much larger than unity as it is as- 
sumed in our model. In this case, high-order corrections 
are expected to be about l/£ or r\ n times, i. e. » 30-40 % 
of, our leading-order result. For the benefit of a quantita- 
tive analysis, the exact expressions (J2j) of the coefficients 
Cj. n have to be used in Eq. (TTJ). 

By employing the laser parameters of ELI (70 PW 
TkSapphire laser focused onto 0.8 /*m, 5 fs pulse dura- 
tion at 1 Hz repetition rate) we would obtain up to two 
photons per pulse from two-photon merging and even the 
merging of ten laser photons would be experimentally ac- 
cessible (two photons per hour). 

In the collision of a laser pulse with a proton beam 
Delbriick scattering can occur as well. While this process 
is clearly distinguishable from the process at hand inso- 
far that the frequency of the scattered photon is different 
from uj n for any n, the question arises if the Delbriick 
scattering can compete with the photon merging process 
in terms of probability. For the considered setup, the an- 
swer is no. In fact, starting from the Born-approximated 
cross section of the Delbriick scattering which is certainly 
suitable at Z = 1 (see, e. g., Q), it can be shown that 
the order of magnitude of the ratio between the number 
of photons scattered via Delbriick scattering and the rate 
TZi is given by (cv/£) 2 . This quantity is very small in the 
situations considered above. 

In the non-perturbative regime of laser-proton inter- 
action when Xn ~ 1 , electron-positron pair creation may 
occur. We show that the depletion of the laser energy due 
to this process is negligible. Thus, in the direct channel 
the pairs are created as a tunneling process when £ 1 
[l9| . The order of magnitude of the pair production rate 
at £ 1 and at laser intensities close to I cr (in the rest 
frame of the proton) is roughly a 2 m. For the data in the 
first (second) example considered, about 10 6 -10 7 (10 8 ) 
pairs are produced per pulse corresponding only to a few 
/iJ (few tens of /zJ) energy loss which is negligible with 
respect to the pulse energy of 3 J (1.8 J). 

The high-energy photons produced via VPEs still 
travel through the laser beam almost in the opposite di- 
rection and could decay via the production of electron- 
positron pairs. The rate of pair production by a photon 
with energy ui in a laser field with £ ^> 1 and xi ~ 1 can 
be estimated as ~ 10 _2 am 2 /wi by employing the results 
in Then, the lifetime in the laser field of a photon 
with an energy ~ 400-500 MeV is around 14-17 fs, i. e. 
longer than the considered laser beam durations. There- 
fore, the photons produced via VPEs are able to escape 
from the region where the strong laser field is present and 
can in principle be detected. 

In conclusion, we have shown that existing proton 
accelerators combined with next generation table-top 
petawatt lasers permit to measure non-perturbative re- 
fractive VPEs in a laser field. These become manifest 



in the opening of multiphoton channels of interaction as 
well as in the enhanced contribution of high-order effects. 
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